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ABSTRACT 



A baroclinic model programmed for numerical solution is used to in- 
vestigate the effects of diabatic heating on the development of pressure 
waves for periods up to 24 hours. Several forms of the model are studied 
with results being obtained for various wavelengths. The results lead to 
some meteorological improbabilities; however, the relative development of 
the pressure wave appears valid for short time periods. 
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Introduction 



Most numerical prediction methods for obtaining prognostic charts of 
the pressure field neglect diabatic heating effects under the assumption 
that their influence is negligible for short period prognosis* 

Under certain circumstances, however, the exchange of sensible heat 
between an air mass and the underlying surface may be considerable. It is 
therefore the purpose of this paper to investigate the effects of such 
diabatic exchange on the development of pressure waves for periods up to 
24 hours using a numerical model* 
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1. 



The Model 



The model with which the aforementioned heating effects were studied 
was developed from the geostrophic relationship, the vorticity equation, 
and the thermodynamic equation, 

V=|-XViZ^; (1) 

'o 

ill + V-V('f'tS)=foiLk:2 , (2) 

^ + V-V^ + 5CJ = - ^ . (3) 

ar t 

The dependent variables are expressed as the sum of a base state plus 
perturbation quantity (indicated by a prime). 



<Z5= 0 


(4) 


y = U Cf)e + v'(y, ' 3 ,T,'fc ); 


(5) 


6J = ^' ( 


(6) 


Q = Q C ?) -t- Q'( p) QCp) =o, 


(7) 




(8) 



Equations (2) and (3) are then linearized by the usual method of assum- 
ing that they are satisfied by the basic quantities alone and by the total 
quantities (sum of the basic plus perturbation). Subtracting the equations 
for the basic flow from the total equations and neglecting products of 
perturbation quantities yields 

, (9) 

i X ax d? 

^ ■+ ta. ^ ■+ O' a 0 + 6 oj* “ ~ -Q (10) 
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Equation (10) is further modified by substituting for the third 

term 

^ ^ ^ (11) 
ax c&f> 

yielding, 

^' + + s w‘ = . (12) 

ifit >X o(f f 

A diagnostic equation for the vertical velocity is obtained by 
differentiating (9) with respect to p, taking the Laplacian of (12), and 
subtracting, resulting in 

3i^ U a Pj" . (13) 

c 4 p ax apJx ap"^ f 

Solutions are now assumed to be of the form 



z ' ° ^ = ACp,'6’)sii^K)<+'B(pj“6‘)ciojkxj 


(14) 






= C C|0,t) 5ii^ )\x +J>Cp,'t) coi'kXj 


(15) 


Q = Q, r s 1 ^ r c Pi 4- ) . 


(16) 



Next equations (14) through (16) are substituted into equations (9) 
and (13), and the coefficients of cos kx equated giving 

^ h ^ C , ) (17) ’ 

--^35^ — t'C p (18) 

- g ^ 3 ,g(a SIK C - O CO^ S. ); (19) 

C p' 

^ 5 AJl ^ +X (20) 

f,'- p" 



3 



Equations (17) through (20) are then solved numerically as an initial 
value problem as outlined by Rosenthal [ij . Given initial values of 
A and B, equations (19) and (20) are solved to obtain C and D. These in 
turn are substituted into (17) and (18) to obtain prognostic values for 
A and B. The computational cycle is continued for the desired time in- 
terval. 
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3. 



Procedures 



Using finite differencing, the diagnostic equations, (19) and (20), 
can be approximated at the odd levels as 

4.) S; - jg Itap^ (^ - 3 ... ) (I9a) 

*’’ «^S,) S - H. ; 

= CUu, -Jf +J-<3 hAp#rA. -A 1 

r V L -/ ‘ 7=.' 51’ ^ I 4 - ‘ 

(-6 

+ -L Cof -I- = - H' . (20a) 



i = 3,5, 19 

At the even levels C and D are obtained by simple averaging, 



c*-S4±£i^ 



1>. = 

L 



SiL+2t;:l . 



(21) 



< = 2,4, 



,20 



Boundary conditions are that C and D at the top and bottom levels 
are zero, implying a vanishing vertical velocity there. 

The prognostic equations, (If) and (JL6) are approximated at the even 
levels by 

liSil " C.., ) 5 G, , (17a) 

h r a A p 

AJx ^ (X. ~J> ) = G-* . ( 18 a) 

c ~ 2,4, 20 

At the odd levels, A and B are obtained by averaging. 






C+l 






( 22 ) 



t= 3,5,, 
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A vertical increment of 50 mb is used, affording considerable vertical 
resolution. This yields a 21 point vertical grid with 0 mb as level 1 and 
1000 mb as level 21, as shown in figure 1. 



Figure 1 - Finite Difference Grid 
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The values of C and D are obtained by a simple relaxation scheme as 
follows; 



m 

n 


6 ^ w- A 


(23) 


n = 


•<'; - £. 


(24) 


• 


♦ --A 




F. = 

L 


f;:!. , 


(25) 


19 , and 


f;* so, then 








( 26 ) 
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(27) 



^: = + f;'- 

6- /«?, 3 

Equations (17a) and (18a) are solved by an iteration technique as 
follows, 

)=A- ('■^)+A£ C(5^ C-t) +G/'^Vt+A-t) ], 

A ^ 

=s^. CG-/ 

A 

C = 2,h,...20, n = 1,2 

Estimates at time t + ^ t are indicated by the superscript n while 
final estimates of quantities at time t are not superscripted. Initial 
estimates for and G* tJ are 



(28) 

(29) 



(Sr/'^ (-ii-A-t) = Q- C-t)j (30) 

^ C-^J. (31) 



The iterative process is continued, recalculating C and D for each 
new value of the superscript, until 



1 A- (-ftzjif)- I ^ 'C i 



infi) 



I ^ -r. 



(32) 

(33) 



When the above conditions are fulfilled, the new prognostic values 
of A and B are then used to calculate corresponding values of C and D. 
The process is continued for the desired time interval. 

The amplitude of the pressure wave and phase angle can then be 
calculated using the relationships, 

z’ = 'Ra ( cos hx - ^ , (34) 
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R * 
z 


-J a- -i- 

V A + B, 


<35) 


tan ^ 


= A/B. 


<36) 



The mod^l was ptogranuned for solution by the CDC 1604 computeir #f the 
U. S* Naval Postgraduate School. The time stepi t^ was taken as 30 
minutes, and the tolerance, > was 1 cm* Solutions for the amplitude 

and phase angle were obtained for every 6 hours up to 48 hours. 

The heating function was introduced in the layers below 500 mb as a 
function of pressure in a manner suggested by Haltiner [ 2 ] . For this 
study, v?as set equal to zero and 

Q->. - ^ (37) 

T S'OO 

In 07), represents the amount of heating being considered in 

joules/gr sac, p is the particular pressure level under consideration 
below 500 mb, and p is the 500-mb pressure. For the purpose of this 

investigation, was taken as 4 x lO"^^ joules/gr sec which is approxi- 

mately equivalent to heating of 400 langlays/day* Haltiner and Wang fsj 
propose that the exchange of heat between an air masa and underlying suirfaee 
in excess of 800 langleys/day« However, it is felt that the 800f 
figure represents an extreme, and that the value of 400 langleys/day is more 
representative^ 

The model is such that the heating function is constant in time and 
space for a given problem, while the phase angle and amplitude ©f the pres- 
sure wave are set initially and then allowed to develop. 

Various permutations of the model were tested with different phase 
angles for the heating function^ and with different wave lengths. 
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In all cases, S was set equal to .5 x 10"^ cgs units, while /3 was 
-13 -1 -1 

equal to 1.6 x 10 cm sec . The wave lengths considered were 2000, 
6000, and 10000 kilometers. 
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Results. 



The results are expressed as a comparison between a no heat case and 
heating introduced with various phase angles. 

a. The first model studied had linearized vertical shear with 

= 40 m/sec and - f ^ , L = 50, 100. ..1000. The initial 

000 

value of B at all levels was 97.4 metersf and A at all levels was zero. 

The phase angle of the pressure wave was therefore zero at all levels. 

With the initial input of heat in phase with the pressure wave, the 
development at the end of a 24-hr period in the lower layers was consider- 
able for all wavelengths, (figure 2). However, in the layer 850-700 mb, a 
decrease in amplitude appears for the very long wave, 10000 km. The long 
wave, 6000 km, showed a similar pattern of development with a minimum of 
development in the same layer, 850-700 mb. 

The movement of the pressure wave, as deduced from the phase angle, 
shows that the heating effect caused an increased movement of the long 
waves in the lower layers, (figure 3). Without heating, the very long 
wave retrogressed slowly at all levels throughout the period. For the 
initial 6-hour period with heating, it retrogressed slowly at all levels, 
then quite rapidly in the lower levels. As seen in figure 3, at the end 
of the 24-hr period, the retrogression in the upper levels is less than 
without heating, while in the lower levels it is considerably more. In 
24 hours, the long wave progressed approximately one-quarter wavelength 
without heating. With heating, it retrogressed initially in the upper 
layers, while progressing in the lower. After 14 hours it progressed at 
all levels, so that after 24 hours, the net effect is progression. The 
short wave progressed over three-quarters of a wavelength without heating 



10 



at all levels, and slightly less with heating. 

Choosing the 500-mb level as representative of the upper levels, a 
gradual increase in amplitude with time and heating was noted for all 
wavelengths, (figure 4). 

The same model with the heating function initially 180 degrees out 
of phase with the pressure wave, showed that at the end of 24 hours, the 
two waves increased in amplitude below 750 mb, while all wavelengths 
increased below 900 mb, (figure 5). The 10000-km wave was the only one, 
when compared to the no heat case, which showed a larger amplitude at all 
levels. The 6000-km wave exhibited a marked decrease at all upper levels, 
while the short wave decreased slightly. 

The very long wave retrogressed in all but the lowest layers. The 
short wave increased its progressive movement to approximately one wave- 
length, while the long wave decreased its movement, progressing slowly in 
the lower layers, and retrogressing in the upper layers, (figure 6). 

A look at figure 7 for the development at 500-mb clarifies the case 
of the very long wave. It too decreased in amplitude until after approxi- 
mately 18 hours had elapsed. This was true for all other levels. 

For initial phase angles of the heating function of 90 and 270 degrees, 
the same marked increase in amplitude in the lower levels is evident, (fig- 
ures 8 and 9). In the upper levels, for an initial angle of 90 degrees, 
the short wave decreased nearly uniformly, the long waves increased. For 
an initial angle of 270 degrees, the opposite appeared. 

b. A second variation of the basic model included a linear 
vertical shear to 500 mb, and then a parabolic variation of the zonal 
wind to the top of the atmosphere. In this case, equaled 20 m/sec, 

and 30 m/sec. A further modification was to allow for tilt of 
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the pressure wave by varying the phase angle in accordance with the 
relationship, phase angle equals q SO looo » 

In all cases, this model showed extreme development below 900 mb for 
all wavelengths, regardless of the initial heating angle. 

When the heating phase angle was zero, the 10000-km wave developed 
very slightly at upper levels, and decreased in the layer 650-800 mb, (fig- 
ure 10). The 6000-km wave developed well above 700 mb, but with a minimum 
of development at 800 mb. The short wave evinced a nearly uniform increase 
in development at all levels above 800 mb. 

For heating with an initial 180-degree phase angle, the 2000-km wave 

showed slight additional development above 700 mb except for a slight de- 

crease in the layer 300-200 mb, and the 10000-km case showed increased 
development throughout, (figure 11). Again, an examination of the 500- 
mb level reveals that the very long wave did decrease in amplitude until 
beyond 12 hours, (figure 12). 

A 270-degree heating phase angle produced increased amplification for 
all levels for the short wave, while both long waves decreased above 800 
mb, (figure 13). 

A 90-degree heating phase angle resulted in an increased amplitude of 
the 6000-km wave at all levels, (figure 14). The 10000-km wave showed a 
decreased amplitude from 900-800 mb, and increased above 800 mb. The 

short wave had decreased development at all levels above 900 mb. 

c. A number of other cases were studied, however, the results 
were comparable to those reported, and therefore are not shown. 
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5. 



Conclusions. 



The restrictions imposed upon the basic model lead t© certain 
meteorological improbabilities, i.e., excessive speed of movement in 
certain cases and over amplification after twenty^four hours. Howeverj 
it is felt that for short time periods, the relative development in the 
layers above 900 mb is plausible even on a quantitative basis. 
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APPENDIX I 



Figures 2 through 14 referred to in text are included here for 
easy reference. 
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